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Abstract 

Quantum dynamics of the electron wave function on one-dimensional lattice is considered. The 
lattice consists of N equal sites and one impurity site. The impurity site differs from other sites 
by the on-site electron energy E and the hopping integral C. The wave function is located on the 
impurity site at t = 0. The wave packet is formed which travels along the lattice, and reflects from 
its end. Reflections happen many times (Loschmidt echo) and this phenomenon is considered in the 
second part of the paper. Analytical expressions for the wave packet front propagation at different 
values E and C are derived, and they are in excellent agreement with the numerical simulation. 
The obtained results can help in interpretation of recent experiments on highly efficient charge 
transport in synthetic olygonucleotides. 
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I. INTRODUCTION 



Charge transfer (CT) is the key event in many processes in hving and inorganic matter. 



The particular interest has 
olygonucleotides was found 



3een expressed after the effective charge transfer in synthetic 
5|. There are possibihties to synthesize the DNA analogues 



with the bases sequences up to dozens base pairs. For example, the double strand consisting 



of 100 base pairs adenosine-thymine (A-T)ioo was synthesized 



6|]. This synthetic olygonu- 



cleotide have regular structure and very effective with regard to the CT. Synthetic DNA 
and polypeptides can be utihzed in nanobiology [?, They are also considered as poten- 
tial molecular wires 
process |lO |. 



9.". It is worth noting that the CT can occur as the one-step coherent 



Experiments on the charge transfer are usually done when donor D and acceptor A are 
attached to both sides of the molecular chain: D — {b)^ — A. The charge is transferred 
from donor on the chain and travels along the chain —{b)^— consisting of N repeating 
DNA bases (usually adenine). The charge can be trapped by an acceptor and registered by 
photophysical response (fluorescence) or electrochemical reactions. 

If an acceptor is absent then the charge is locked on the chain and can reflect from 
the end. Many reflections can happen (the phenomenon of Loschmidt echo, i.e. multiple 
returning to the initial state, was extensively studied in classical and quantum systems 
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15|). V. Benderskii with colleagues investigated the analogous quantum dynamics of 



vibrational excitations in ID molecular chain 
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18| . It worth noting that the property to 



form the moving wave packet is an attribute of discrete lattices. The diffusive spreading of 
the initially localized excitation is usually observed in continuous models. 

The recurrence phenomena to the initial state are well known. This property is also typical 
for some dynamical systems. The Fermi-Pasta-Ulam recurrence is an example 
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20|. 



The other example is the interaction of the isolated vibrational mode with the continuous 
spectrum |21 |. 

The recurrence in the quantum systems was for the first time considered by R. Zwanzig 



22| for the discrete equidistante spectrum when single level is initially populated. The 



analysis of the Zwanzig' s problem is done in 16|. The same authors generalized their 
consideration in the next paper jl^]. The quantum dynamics of the vibrational excitation 



on the one-dimensional lattice was considered in 



y. 



It was shown that if the excitation 
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is initially localized in the lattice center, then the energy is not distributed homogeneously 
along the lattice but conversely it localizes and reflects from the lattice ends many times. 

In the present paper we consider an analogous problem of the wave packet propagation on 
the lattice and multiple reflections from the lattice ends. One of the goals is an explanation 
of efficient charge transport in synthetic olygonucleotides. 



II. SETTING UP A PROBLEM 



We consider a lattice consisting of N identical sites and one impurity site at the left 
lattice end. The excitation is initially confined to the impurity site. An excitation can be 
electronic or vibronic. The electronic excitation is considered for definiteness. 

The wave function amplitude on the impurity site is labelled by a{t), and amplitudes 
on other sites bi{t) {i = 1,2, .. . ,N). The electron can hop on the neighboring site. The 
electron energy on all sites except the impurity site is zero, what corresponds to the choice 
of the reference point for electron energy. The hopping integral C — 1, what corresponds to 
the choice of energy unit. E is the on-site energy of the impurity site, and C - the hopping 
integral between the impurity site and the nearest site of the lattice. The challenge is to 
find how the electronic populations of all sites change in time. 

The hamiltonian in the {N + 1) x [N + 1) matrix representation reads: 
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(1) 



with the wave function 



^{t)^ait),bi{t),b2it),...,bNit). 
In the second quantization representation the hamiltonian is 



(2) 



Af-l 



N 



H = Ea^a- + J] 6+ ,6r + ^ bt^br + C {a+b^ + fe+a") 



(3) 
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FIG. 1: The dependence of the wave function amphtude a{t) on the impurity site vs. time. 
The wave function returns many times (four times in the figure) to the impurity site with shghtly 
decreasing amplitudes. Parameters: E = 0, C = \/2, N = 50. 



The dimensionless {h = 1) Schrodinger equation for the wave function has the form 

ia = Ea + Chi 
ibi = 62 + C'o 

4=61 + 63 (4) 



ib 



N 



JN-1 



with the initial conditions a{t)\^^Q = 1 and bi{t)\^^Q = (i = 1, 2, . . . , A^). The formulated 
problem mimics the experiments on the charge transfer in synthetic DNA {2]. 

If equations (jl]) are integrated numerically then the phenomenon of repeated reflections 
from the lattice ends is observed, i.e. Loschmidt echo or the electronic "ping-pong" (the 



23|). The result is shown in Fig. [T] 



electronic ping-pong was found experimentally 

The goal of the present paper is the consideration of the wave packet propagation and 
its first reflection from the lattice end. This time interval corresponds to the time range 
< t < 100 in Fig. [TJ Next part of the paper considers the multiple reflections of the wave 
packet. 

The quantum dynamical problem of the wave packet evolution can be treated as the 



interaction of the impurity site with the reservoir. And it is solved using the expansion in 
terms of the eigenfunctions of the reservoir. 

The lattice without the impurity site (reservoir) is described by the NxN tridiagonal ma- 
trix with zero leading diagonal and the unity values on the secondary diagonals. Eigenvalues 
e{k) and eigenfunctions qi{k) of this matrix are well known: 

Tik 



e{k) = 2 cos 



sm 



vr/c 



(5) 



It is more convenient to consider the problem in terms of amplitudes of modes h{k) instead 
of amplitudes on the lattice sites 6j. These values are related by the orthogonal relationship: 

N N 

b{k,t) = Y,^r{k)h{t)- hit) = Y,qiik)b{k,t). (6) 



i=l 



k=l 



Then the system can be written in the form: 



TV 



ia{t) = Ea{t) + C ^ qi{k) b{k, t) 



(7) 



k=l 

ib{k, t) = e{k) b{k) + Ca{t) qi{k) 
Amplitudes b{k,t) in are expressed through the amphtude a(t): 

b{k,t) = qi{k) exp[—ie{k)t] / a{T) exp[i e{k) rjdr . (8) 

Jo 

Substituting this expression into equation (I7j) for a(t), one can get the integro-differential 
equation for the amplitude a{t) on the impurity site: 



d{t) = -t E a{t) - [ (t - t') a{t')dt' , 

Jo 

where the kernel is given by the following sum: 



N 



1^ 

k=l 



sm 



Tck 



N + 1 



exp 



—2i cos 



TT 



k 



N+1 



(9) 



(10) 



Our primary goal is the solution of ([9]) for the amplitude a{t) on the impurity site. 



III. THE DECAY OF THE INITIAL STATE a{t = 0) = 1 IN THE SEMI-INFINITE 
LATTICE 

Initially we consider the decay kinetics of the initial state a(t)|(=o = 1 iii the semi-infinite 
lattice (A^ — )■ oo). The amplitude on the impurity site is labelled by ao(t) in this case. In 
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the expression ( ITOj) for the kernel B^{t) we use the hmit N ^ oc. The hmiting value (at 
—7- oo) of the expression for B^{t) is labelled by Bo{t). 
The sum ( llUp in the limit — )■ oo is modified to the integral 24| : 

Bo{t) = \im B^{t) 1^^^ = Jo(2t) + J2(2t). 

And in this limit we have the following equation for the amplitude ao(t): 

ao{t) = -lEaoit) - C'^ f Bo{t - t') ao{t') dt'. 

Jo 



(11) 



(12) 



This equation can be solved using the Laplace transformation: 

""'^^^ = p + zE + C'Bo{p) ' 
where the Laplace transform Bo{p) of the function -Bo(^) t24| : 



p 



(13) 



(14) 



The amplitude ao(t) can be obtained by the inverse Laplace transformation: 

d+ioo 



ao{t) 



exp{pt) 



d—ioo 



p + iE + C^^ (v/p' + 4 



p 



dp, d> 0. 



(15) 



The obtained form of the expression for ao(t) is very inconvenient for the analysis and 
numerical computations (integral with the upper infinite limit converges very slowly). The 
convenient form this integral obtains if the integration contour is closed around the cut 



[— 2i,+2i], and the square root a/p^ + 4 should be rewritten in the form: 



—i The integration contour can be closed if amplitude ao{p) has no poles. 

Assuming that poles are absent, close the integration contour and summing the integrals 
along both banks of the cut, one gets the following expression for amplitude ao(t): 



o-oi^) = — / du exp{iut) Im u + E (u + i\/4 — 

J -2 L 2 V 



(16) 



Few special cases can be noted, when the integral (116 
amplitude ao(t) is expressed through the Bessel functions 



is calculated explicitly and the 



ao{t) =Jo(2t); 

ao(t) =M2t) + J2(2t); 

ao{t) =Jo(2t) -zJi(2t); 



{E = 0, C = V2) 
{E = 0, C=l) 
iE = l, C = l) 



(17) 
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FIG. 2: The dependence of the amphtude ao(t), calculated according to vs. time at ii^ = 0. 
Solid Hne: = 0.1; dotted line: = 0.5; dash-dotted line: = 1; dashed line: = 2. Empty 
circles - accurate result for = 2, when ao(t) = Jo{2t). 

An integro-differential equation for the amplitude ao{t) is derived in Appendix A, and 
its solution is obtained as a series in terms of the Bessel functions. These series have simple 
form in few particular cases when E = or C = 1. 

Fig. [2] shows the dependence of the amplitude ao{t) vs. time at £" = and at different 
values of the hopping integral C^. The comparison with the accurate result is also shown 
("accurate" result = numerical integration of equations @)- 

The dependence of the amplitude vs. time in (fT6l) can be estimated by the Fermi's golden 
rule, when the perturbation approximation by the small parameter is used, and when an 
exponential decay is valid: 



More accurate expression for the amplitude ao{t) at small values of parameter C is given in 
Appendix B. 

Now we consider the case when amplitude ao(p) has poles. An analysis shows that one 
pole exists in the region E > 2 — C^, and two poles - when E < C"^ — 2. There are no poles 
in the region \E\ < 2 — (C^ < 2). Fig. [3] shows number of poles vs. values of parameters 




(18) 



E and C. 
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FIG. 3: Number of poles of the integrand of function (|16p in regions divided by dashed hnes. 



As an example we consider the case of one pole in the region E > 2 — C"^ (C^ < 2). The 
pole is located at the point p = —is {e > 2). The relation between parameters E and C 
follows from ( |T5l) : 



The overall amplitude is the sum of the main term (fT6|) and the pole contribution (!20|) . 
The dependence of amplitude vs. time when one pole exists, is shown in Fig. HI In this case 
the amplitude on the impurity site does not decrease to zero and ao(^)|f_^oo ~ 0.75. 

In Appendix C it is shown that ii \E\ > 2 — (C^ < 2), then there exists the bounded 
state localized at the lattice end. The value e (see f|T9|) ) is the energy of this localized state 
and its contribution to the amplitude a^lt) coincides with the contribution from the pole 
summand. Note also that the integral term f|T6l) is the contribution to the amplitude from 
the continuous spectrum which lies in the range [—2, 2]. The fraction of the initial state 
is trapped by this localized state and the smaller part of the initial state goes out to the 
reservoir. The effect of returning to the initial state also decreases. This case (existence of 
a pole) seams to be less interesting and we consider only an absence of the localized state 
when the integral (|T6|l is the accurate result for the amplitude ao(t). 




(19) 



The contribution from the pole Aoq is equal to 



exp{—iet) 



(20) 
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FIG. 4: The dependence of amplitude ao{t) vs. time in the presence of one pole. Dashed line - 
contribution to amplitude from the integral summand (jl6p : dots - modulus of the pole summand 
(equal to 0.75) (j20p : empty circles ~ numerical result. Parameters: E = 2, C = 1 {e = 2.5). 



Such a detailed analysis of the amplitude ao(t) decay in the infinite lattice is done by 
two reasons. First, the decay kinetics in the infinite lattice coincides with the kinetic in the 
finite lattice consisting of sites in the time range < t ^ N. Second, as will be shown in 
the second part of the paper, the overall amplitude in the finite lattice can be represented 
as a sum of partial amplitudes, first of which is just the amplitude ao{t). 

IV. WAVE PACKET PROPAGATION ON THE LATTICE 

The dynamics of the wave function on the lattice, i.e. spatiotemporal evolution of ampli- 
tudes bj{t), is considered in this section. Initially the expected spreading of the wave packet 
is observed. But then amplitudes on the sites with larger numbers increase. And when the 
amplitude a{t) becomes small, the wave function forms the well defined wave packet with 
the sharp forward front. This behavior is shown in Fig. O 

An expression for amplitudes bj(t) is obtained by the substitution expression ([8]) for the 
mode amplitudes b{k,t) into expression (Q. Then one gets the following equation 

= ~WTl ^ {wTl) {wTl-^^ exp[-te{k)t] exp[-i£(A;)r]a(r)dr. (21) 
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FIG. 5: Spatiotemporal evolution of the wave packet. Solid line: t = 2; dashed line: t = 5; dotted 
line: t = 10; dash-dotted line: t = 15. Site j = corresponds to |ao(t)p. Parameters: E = 1, 
C2 = 0.5. N = 100. 



Initially we consider the case of semi-infinite lattice (A^ — t- oo). Changing the summation 
in (pTi) by integration and replacing a{t) by ao(t), the following expression can be obtained: 

i2c n /"* 

^i(^) ~ / dksm{k) sm{kj) exp[—ie{k)t] / exp[— i£(A;)r]a(r)dr , (22) 

Jo Jo 

where e{k) = 2cos(A;). 

Consider the case when time t is large and amplitude ao{t) is very small, |ao(t)| ^ 1. 
Then the upper limit in integration over time in (1221) is infinity. The obtained integral is 
the Laplace transformation ao{p) taken at p = —ie{k). As a result we have: 

Expression fl23l) can formally be generalized on the negative values of index j. Aiming 
this in mind, the multiplier sin{kj) is represented as the difference of two exponents. Then 
expression f l2^ describes the superposition of two wave packets freely travelling to the left 
and to the right in the infinite (to both sides) lattice. Every wave packet is normalized to 
unity. Consider separately the wave packet propagating to the right. Label this wave packet 
by b'f: 

LiMu\ J, sin(fc)exp{-z[A;j + 2cos(A;)t]} 



hfit) = / dk ' \ ^ ' -^^^ tHt^' -oo < j < +00 . (24) 

^ 71 Jo E - 2 cos{k) + exp{-ik) ^ ^ 
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FIG. 6: Snapshots of the impulse in semi- infinite lattice at t = 25 (solid line); t = 50 (dotted 
line); t = 75 (dashed line). Calculations according to (j24p . Empty circles - numerical integration 
of ([4]). The coincidence of the numerical and analytical results is also excellent for t = 50 and 
t = 75 (mean square error (MSE) < 10~^). 

The considered impulse bj (in the region j > 0) is the following difference: 

6,(t)=6f(t)-rj(t). (25) 

The addend 6™j(t) is the "tail" of the impulse bj{t), and it is small at large times, e.g. ~ 
Then it can be assumed that U^^{t) is a good approximation for hj{t). The snapshots of the 
impulse at different time points are shown in Fig. [61 One can see that the constrain by the 
impulse 6™^(t) is a very good approximation. 

Expression &™^(t) (see ( 124]) ) is a reliable approximation for the propagating impulse on 
the time interval t < N/2 for the finite but comparatively long lattices. Moreover, as will 
be demonstrated below, this expression allows to describe the impulse reflection from the 
lattice end on the time interval N/2 < N. 

Lets analyze the expression (p^ at large times. As the integrand is the fast oscillating 
function, the stationary phase method can be applied. If j < t/2, then two stationary points 
exist which are determined by the equality sin(A;) = j/2t. Then the following expression can 
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FIG. 7: The impulse front at t = 200 {E = I, C"^ = 2, N = 500). Dotted line - numerical result, 
solid line - the asymptotic expansion (j28]) . 



be obtained: 



-iC 



^y Tct cos{ko) 



[/(tt — ko) exp(i7r/4) + /(tt — k) exp(— z7r/4)] , ko = arcsin(j/2t). 



(26) 



Here f{k) is the integrand in expression (j2l]). 

Fig. [7] compares the expression fl26|) . obtained by the stationary phase method, with the 
accurate answer. 

If j ~ 2t, then both stationary points lie close to the point k = 7i/2 and the stationary 
phase method gives false results (if both points coincide then the answer diverges). Lets 
consider this region {j ^ 2t) in more details. In this case the expression for amplitudes ([2] 
can be simplified: 

1 



iC exp(— O.Sttj) 



exp 



dx- 



j-2t X 3 
3 



. ,27) 



7r(2 - C2) z-x ' ' ' 2 - C2 

If \z{t)\ ^ 1 (what is valid for t ^ 1), then the asymptotic expansion of ( l26l) can be 
rewritten as a series in Airy functions: 



6,(t) 



2iC 



zm-c^) 



Ai(j) 



^ -Ai'(j)-4TTAi"(j) + --- 



3 



J-'2t 

rt ■ 



(28) 



iz{t) z^{ty 

Thereafter an approximate expression for the impulse front can be written. Using the 
asymptotics of the Airy function, one can notice that the dominant term, describing the 
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FIG. 8: The impulse front at t = 200 {E = 1,C^ = 2). Empty circles - numerical result, solid 
line - asymptotic expansion (j28p . MSE < 10~^. 

impulse front, has the form: 



Fig. |8] shows the comparison of the asymptotic expansion ( l28l) with the numerical answer. 
Expression ( l29l) allows to describe qualitatively the form of the impulse front. The depen- 
dence on the lattice number j is defined by the argument (j — 2t)/\/t (see ( 1271) ). From ( 128|1 



decreases oc 1/ \/t. The front velocity is 2 - maximal possible group velocity. The front has 
the sharp profile. Thus, if the finite lattice is considered, then until the front achieves the 
lattice end, its propagation is the same as in infinite lattice. Front of the impulse achieves 
the lattice end at t ~ N/2. After that the impulse refiects from the lattice end and the front 
moves back to the impurity site. 

Now we consider the impulse refiection. Lets assume that the lattice is long enough. 
Then amplitude a(t) is negligible and impulse is totally formed. The impulse front moving 
with the velocity f = 2, reaches the lattice end at t ~ N/2. After refiection it moves with 
the same velocity f = 2 in the opposite direction. And at t ~ the impulse reaches the 
impurity site on the left end. Amplitude a{t) starts to increase. We consider the impulse 
evolution at time range t < N, when amplitude a{t) is small and when it coincides with the 




(29) 



one can see that the impulse front slowly spreads: its width increases oc \/t and amplitude 
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amplitude ao{t) — amplitude in the infinite lattice. 



V. IMPULSE REFLECTION 

To derive an expression for the reflected impulse we come back to the expression ( |2T]) for 
amplitudes bj{t). We assume that the time is large such that amplitude a{t) is small. Then 
the upper limit in integration over time is infinity and the substitution a{t) — )■ ao(t) is made. 
As earlier, the obtained integral can be represented as the Laplace transformation of ao(p) 
at p = —ie{k). Then the following expression can be deduced: 

N 



hj{t) = -— — ^ ^sin(A;) sin(A;j) exp[-22t cos(A;)] (E - 2cos(A;) + C'^ exp{-ik)j , (30) 
~'~ k=i 

where k = T[k/{N + 1). 

This expression can be modified using the Poisson summation formula: 

7 OO 

K=oo m,=oo „ 

f{k) = / f{x) exp{—2Tcimx) dx (31) 

fc=— oo m=— oo_^ 

and the function f{k) differs from on the interval [0, 1]. After substitution of variables 
TTx/ {N + 1) = y, the following series can be obtained: 

_ 2C f smysm{jy)exp{-i[2tcosy + 2my{N + 1)]} 

hi^)- vr ^ 7 E -2cos y + C^exp{-zy) "^^^ ^^^^ 

This series can be represented as a sum of two series, expressed through the amplitude 6"^ 
— the amplitude in the infinite lattice at large time (see (12^ ): 

m=oo m=oo 

hit)= E ^lt2miN+l)it) - E ^-1+2„.(iV+l)(^) (33) 
m=— oo m=— oo 

For the separating terms, which are significant in the considered time range (1 ^ t < A^), 
the expression flMl) should be analyzed and the terms, where the phase has the stationary 
point, should be found. From the first sum the single term with m = is left which describes 
the incident impulse 6™^ From the second sum - the term with m = 1, i.e. &™j^2(iv+i) 
describing the reflected impulse. Thus, if t < A^ then the impulse can be represented as the 
sum of two terms: 

6,(t) = 6f(t)-6!_i+2(^+i)(t). (34) 
14 
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FIG. 9: Reflection of the impulse from the lattice end (A^ = 100). Time t = 75. Parameters: 
E = 1, = 0.5. Computation according to (|26l) . Solid line - reflected impulse, dotted line - tail 
of the incident impulse. At i = 75 the impulse front, moving with v = 2, passed ~ 150 lattice sites: 
100 sites before reflection and 50 sites after reflection. An excellent agreement with numerical 
simulation is observed (MSE < 10~^, numerical data are not shown). 

Fig. [9] shows the impulse reflection calculated according to ( 134|) . 

At t ~ the impulse front returns to the lattice beginning and starts to interact with 
the impurity site. Amplitude a{t) increases. Our approximation concerning the smallness of 
a{t) becomes invalid. The system returns to the initial state. This process will be analyzed 
in the next paper. 

VI. CONCLUSIONS 

The electronic impulse propagates starting from the impurity site at the lattice end. If 
the electronic excitation is entirely located on the impurity site then, depending on the 
impurity site parameters, few scenarios of the wave function dynamics are realized. One, 
most common case, is the formation of impulse and its propagation along the lattice. The 
impulse moves with the maximal group velocity v = 2 and has the steep forward front. After 
reaching the opposite lattice end, the impulse reflects and moves in the opposite direction. 
This process repeats many times. The other scenario is realized, if the parameters of the 
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impurity site are such that the bounded state is formed. Then the wave function is partially 
trapped by this bounded state. 

The analytical expressions describing the impulse dynamics coincide with the accurate 
numerical simulation. 

The problem of the electronic excitation travelling along the lattice, being formulated 
in a rather simple form, can be relevant to recent experiments on the charge transfer in 
synthetic olygonucleotides and polypeptides. 

At one time it seemed that the polaron mechanism can adequately explain the charge 
transport in biomacromolecules. But there exists at least to weak points in this approach. 
First, - polaron formation. After a charge is transferred to the chain it should be more or less 
localized for some time to allow the polaron formation. Otherwise the wave function spreads 
very fast along the lattice because of very small hopping time. Second problem, - how the 
polaron can get an initial momentum necessary for the coherent movement. The discussed 
mechanism of quantum dynamical charge transfer allows to exclude this difficulties in ex- 
plaining the high efficient CT over long distances: the electronic wave packet, responsible 
for the charge transport, spontaneously forms and moves. 



Appendix A 

An expression for amplitude (see (ITB]) ) can be rewritten in the following form: 

-2 

This expression can be rearranged to the differential equation by time. It can be done if the 
denominator of the integrand is represented as the differential operator with the replacement 
u — 7- —id/dt. Then the action of this differential operator on both sides of equation flAip 
eliminates the denominator. The resulting integral is expressed through the Bessel functions. 
As a result the following differential equation is obtained: 

(1 - C^) do + iE{2 - C^) do - (C^ + E^) ao = [Jo(2t) + J2(2t)] . (A2) 

Hence the inhomogeneous linear differential equation of the second order is derived. Initial 
conditions for this equation are: ao(t = 0) = 1 and daoldt\t=o = —iE (the second condition 
follows from equation (fl2|) ). 
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The solution of this equation can be formally written as the integral over time. But 
this integral is difficult to evaluate analytically and numerical calculation is more complex 
compared to the original expression ( TT^ . 

Nevertheless we demonstrate that the solution can be represented through the Bessel 
functions. And there exist few particular cases when the solution can be written in the 
rather compact form. 

If C = 1 then flA2|) is the first order equation. Its solution can be expanded into a series 
in Bessel functions: 

oo 

ao{t) = Jo(2t) - iEJi{2t) + (1 - E-^) ^{-lEf J,,(2t). (A3) 

k=2 

In the case when = the equation (1A2|) has no first derivative and its solution is a 
series in even Bessel functions: 

oo 

ao(t) = Jo(2t) + (2 - 5^(1 - C^ f^' M2t). (A4) 

fc=i 

In the general case the solution of flA2p can be written as the following series: 

k=2 

where 



ao(t) = Jo(2t) - zEM2t) - ^^^^ _ J^ihxl - hx^^) M2t) , (A5) 



'E^ + C 



a 



2-C 



2 



[2(1-C2)_ 



2 



.2-C' 
ai=a- I— — ; a2 = -a^ 

2(1 - C ) (A6) 

4 + a2 

h= ' — ; ^ = 1,2 
+ 2 - aiXi 



X, =i ( ai± + ) ; i = l,2 



Sign for Xj in (]A6P is chosen so that < 1. 



Appendix B 



When the hopping integral C is small then amplitude ao(t) slowly exponentially decays. 
Below we get two approximate expressions for the amplitude when C' <^1. 
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FIG. 10: The transformation of the initial integral (IBlh (dashed line) into two integrals Ii C 
[—2; —2 + ioo], I2 C [+2; +2 + ico\ and the residue in the pole uj = uji. 



The first method consists in modification the integral formula (fT6l) for ao(t). This expres- 
sion can be rewritten as: 

2 

^2 

ao(t) 



27r 



dcj- 



-\/4 — exp(2(X'i(:) 



a; + E)(w + E-a;C2) + C4 



(Bl) 



If C is small then an approximate expression for ao(t) can be written. To do this, the 
integration contour in (IBip should be deformed as shown in Fig. [TOl In the upper semiplane 
the integrand in fIBip has a pole at the point oj = oji. 

1 



(-E{2 -C^) + i C^^^ - AC^ - . 



(B2) 



2(1 -C2) 

When C ^ 1 then the pole lies very close to the real axis (on the distance ~ C^). 
Therefor the contribution of the pole term decreases with time very slowly, - slower then 
the contribution from integrals Ji and I2 (see Fig. ITOll . 

The reason is that integrals Ji and I2 contain the exponentially decaying multiplier 
exp(— cut). Therefor at large time the contribution from these integrals is small in com- 
parison with the contribution from the pole. Thus, at small and taking into account only 
the pole contribution, we get the following approximate expression for the amplitude decay: 



ao(t) ~ exp(za;t) 



4-4C2-£;2- 



(B3) 
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At small the expression for the amplitude can be additionally simplified and the result 
coincides with the Fermi's golden rule: 



ao(t) ~ exp 



-iE - -CV4 - ^2 ) t 



(B4) 



Other approximate expression for the amplitude ao(t) (at ^ 1) can be derived using 
the perturbation theory for the initial equation ( fT2l) : 



do(t) = -iEao{t) -C^ j Boit- t') ao{t') dt'; Bo{t) = Jo(2t) + J2(2t) . (B5) 



Lets introduce the function Hq = aQexp{iEt). For this function we have the following 
equation: 

_ t 
dao 



dt 



-C^ j Bo{t - t') aoit') dt' 



(B6) 



Function oq varies slowly at small C, and in the first order of the perturbation theory we 
can put ao(t') ~ ao(^) in the integrand. As a result the linear differential equation can be 
obtained. Its solution is: 

^ t 



ao(t) = exp 



J{t- t') exp(2Et')[Jo(2t') + J2(2t')] dt' 



(B7) 



Note that this expression is valid when time t is small. If t ^ 1 then expression (1B7P can 
be additionally simplified. Evaluating the corresponding integrals and returning back to the 
function ao{t), we get: 

iEC^ 



ao{t) = exp 



CW4: - E^ 



t-iE{l + ^\ t + 



2V4-E2_ • 



(B8) 



Note that the error estimation of the exponent shows that this error is of the order not 
exceeding t"^/^. 

It also should be pointed out that expression fIBSp for |aoP gives an answer coinciding 
with the Fermi's golden rule (but with extra multiplier exp(C^) in (IBSp ). Fig. [TT] shows 
the comparison of numerical answer for the amplitude ao(t) decay with the approximate 
expressions (IB3P and (IBSp . 
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Appendix C 

Consider now the case when E > 2 — C^, (C^ < 2) and when there exists one locahzed 
bounded state. We find the eigenstate \l/ioc locahzed at the lattice end and solve Schrodinger 
equation -ff^^ioc = ^^E^ioc (for the hamiltonian H see ([I])). The solution has the form of 
exponentially decaying function of j: 

^,oc = Al^,a^-'\, A=(^ + -^) ' . (CI) 



Parameter a is expressed through the energy e: 

a = \{e~^7^^. (C2) 

Energy e is related to parameters E and C by the relationship, coinciding with given by 
(dZD: 



E = e (^1-— J+— v/i^. (C3) 
Expand now the total wave function in terms of the eigenf unctions of the hamiltonian: 

= D "^loc exp(-2et) + C.C.S., (C4) 
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where CCS. means the contribution from the continuous spectrum. The first term is just 
the contribution to amphtude from the locahzed state. It is necessary to find the constant D. 
To do this we consider the scalar product of the locahzed and the complete wave functions. 
Because of the orthogonality of wave functions one can get: ^^loc^^ = -D exp(— iet). In 
this equality we put t = 0. As the wave function is fully localized on the impurity site at 
t = 0, i.e. ^(t = 0) = (1, 0, 0, . . .), then it follows that D = = 1). Thus, the necessary 

contribution to the amplitude on the impurity site from the localized state is: 

a'°' = \%oc{l)fexp{-tet). (C5) 

Substituting the expression \E'ioc(l) = AC^^ from (ICip and expressing A through e (see (IC2p ) 
we get the final answer, coinciding with ( !T7|1 : 

a^t) - ^rf-:'^ , ■ (C6) 
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